- Author:
- aram148 <42922407+aram148@users.noreply.github.com>
- Date:
- 2022-07-19 16:08:46+12:00
- Desc:
- Added CO2 exchange and pH bg models
- Permanent Source URI:
- https://models.cellml.org/workspace/8c0/rawfile/afeeca156d18d9435303b544ea071f7cf868595d/WagnerCano SS Respiration-driven 2020-01-19 - Copy.cellml
<?xml version='1.0' encoding='UTF-8'?>
<model name="my_model" xmlns="http://www.cellml.org/cellml/1.1#" xmlns:cellml="http://www.cellml.org/cellml/1.1#">
<!-- Wagner-Cano SS respiratory model -->
<units name="dim">
<unit units="dimensionless"/>
</units>
<units name="per_s">
<unit exponent="-1" units="second"/>
</units>
<units name="per_L">
<unit exponent="-1" units="liter"/>
</units>
<units name="per_mmol">
<unit exponent="-1" prefix="milli" units="mole"/>
</units>
<units name="per_J">
<unit exponent="-1" units="joule"/>
</units>
<units name="L_per_s">
<unit units="litre"/>
<unit units="per_s"/>
</units>
<units name="L_per_mmol">
<unit units="litre"/>
<unit units="per_mmol"/>
</units>
<units name="mmol_per_s">
<unit prefix="milli" units="mole"/>
<unit units="per_s"/>
</units>
<units name="mmol_per_L">
<unit prefix="milli" units="mole"/>
<unit units="per_L"/>
</units>
<units name="kPa">
<unit units="joule"/>
<unit units="per_L"/>
</units>
<units name="mmol_per_s_per_kPa">
<unit prefix="milli" units="mole"/>
<unit units="per_s"/>
<unit exponent="-1" units="kPa"/>
</units>
<component name="main">
<variable initial_value="0" name="t" units="dim"/>
<variable initial_value="0" name="V_In_air" units="L_per_s"/>
<variable initial_value="1" name="V_In_air_max" units="L_per_s"/>
<variable initial_value="0.2093" name="F_In_o2" units="dim"/>
<variable initial_value="22.4e-3" name="f_o2" units="L_per_mmol"/>
<variable initial_value="101.33" name="P_B" units="kPa"/>
<variable initial_value="0.383" name="V_b" units="L_per_s"/>
<variable initial_value="28.46" name="D_l_o2" units="mmol_per_s_per_kPa"/>
<variable initial_value="0.5692" name="D_m_o2" units="mmol_per_s_per_kPa"/>
<variable initial_value="2.842" name="Q_m_o2max" units="mmol_per_s"/>
<variable initial_value="0.03333" name="p_m_50" units="kPa"/>
<variable name="p_l_o2" units="kPa"/>
<variable name="Q_l_o2" units="mmol_per_s"/>
<variable name="Q_m_o2" units="mmol_per_s"/>
<variable name="p_lc_o2" units="kPa"/>
<variable name="p_mc_o2" units="kPa"/>
<variable name="p_m_o2" units="kPa"/>
<variable name="C_lc_o2" units="mmol_per_L"/>
<variable name="C_mc_o2" units="mmol_per_L"/>
<variable initial_value="8.615" name="C_o2_max" units="mmol_per_L"/>
<variable initial_value="3.573" name="p_50_o2" units="kPa"/>
<variable initial_value="3" name="n" units="dim"/>
<variable name="Ventilation" units="mmol_per_s"/>
<variable name="VentilationPerfusionRatio" units="dimensionless"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>V_In_air</ci>
</apply>
<apply>
<times/>
<cn cellml:units="dim">1.0</cn>
<ci>V_In_air_max</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>Q_l_o2</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>V_In_air</ci>
<apply>
<minus/>
<ci>F_In_o2</ci>
<apply>
<divide/>
<ci>p_l_o2</ci>
<ci>P_B</ci>
</apply>
</apply>
</apply>
<ci>f_o2</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>Q_l_o2</ci>
<apply>
<times/>
<ci>D_l_o2</ci>
<apply>
<minus/>
<ci>p_l_o2</ci>
<ci>p_lc_o2</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Q_l_o2</ci>
<apply>
<times/>
<ci>V_b</ci>
<apply>
<minus/>
<ci>C_lc_o2</ci>
<ci>C_mc_o2</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Q_l_o2</ci>
<ci>Q_m_o2</ci>
</apply>
<apply>
<eq/>
<ci>Q_m_o2</ci>
<apply>
<times/>
<ci>D_m_o2</ci>
<apply>
<minus/>
<ci>p_mc_o2</ci>
<ci>p_m_o2</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>C_lc_o2</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>C_o2_max</ci>
<apply>
<power/>
<ci>p_lc_o2</ci>
<ci>n</ci>
</apply>
</apply>
<apply>
<plus/>
<apply>
<power/>
<ci>p_50_o2</ci>
<ci>n</ci>
</apply>
<apply>
<power/>
<ci>p_lc_o2</ci>
<ci>n</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>C_mc_o2</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>C_o2_max</ci>
<apply>
<power/>
<ci>p_mc_o2</ci>
<ci>n</ci>
</apply>
</apply>
<apply>
<plus/>
<apply>
<power/>
<ci>p_50_o2</ci>
<ci>n</ci>
</apply>
<apply>
<power/>
<ci>p_mc_o2</ci>
<ci>n</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Q_m_o2</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>Q_m_o2max</ci>
<ci>p_m_o2</ci>
</apply>
<apply>
<plus/>
<ci>p_m_50</ci>
<ci>p_m_o2</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Ventilation</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>V_In_air</ci>
<ci>F_In_o2</ci>
</apply>
<ci>f_o2</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>VentilationPerfusionRatio</ci>
<apply>
<divide/>
<ci>Ventilation</ci>
<ci>Q_l_o2</ci>
</apply>
</apply>
</math>
<!--C_lc_o2 = C_o2_max*p_lc_o2/(p_50_o2+p_lc_o2);
C_mc_o2 = C_o2_max*p_mc_o2/(p_50_o2+p_mc_o2); -->
</component>
</model>
